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S ' Abstract 

o, 

O , We present a phenomenological model for the vortex dynamics in 

the peak effect region of weakly pinned superconductors. Wc explain 
the history dependent dynamic response of the metastable vortex states 
^ ' subjected to a transport current and the hysteretic voltage-current 

\l , characteristics observed in the vicinity of peak effect in weakly pinned 

superconductors. A strong variation in voltage current characteristics 

qC ', with the current sweep rate and the anomalous dependence of critical 

f^ ' current density Jc on the magnetic field sweep rate have also been 

^^ I accounted for by this model. 

O 
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In weakly pinned superconductors, the competition between intervortex in- 
teraction and quenched disorder produces a sharp peak in the critical current 
density Jc just below the normal state boundary. This phenomenon, known 
'^ I as the peak effect ||l| occurs when the vortex lattice passes from a state with 

a high degree of spatial order to an amorphous pinned state P, ^. The 

O . peak effect is often accompanied by a marked history dependence in Jc orig- 

inating from different metastable configurations in which the vortex lattice 
couldexist §|,|,|,|. 

Another aspect of the metastability associated with the peak effect is 

5^ I seen in the history dependent dynamic response of the vortex lattice to a 

transport current^]. It is presumed that this behavior is due to the rear- 
rangement of the vortex lattice from one metastable configuration to another 



aided by an external driving force |^]. Moreover, the process of rearrange- 
ment can be halted or in other words, a particular nietastable configuration 
can be frozen in time, by merely switching off the transport currentp]. 
Switching on the current again after a time interval revives the further evo- 
lution of the metastable vortex configuration from where it was halted (long 
term "memory" effect). Further, the evolution of a metastable vortex config- 
uration can be considerably slowed down by slightly lowering the transport 
current!^]. 

Additionally, Zhukov et a/[|lO| have recently reported a magnetic field 
sweep rate (dB/dt) dependent magnetization study in the peak region of 
2H — NbSe2 and in the fishtail region of YBa2Cu^0T. Their results imply 
that in the increasing field cycle dJc/d\dB/dt\ < at a fixed applied field. 
This so called "negative dynamic creep" phenomenon is in contast to the be- 
havior expected from the thermally activated fiux creep where Jc is expected 



to increase logarithmically with the sweep rate {dJc/d\dB/dt\ > 0)|11]. 
However, in the decreasing field cycle, it is found that dJc/d\dB/dt\ > 0. 

In transport experiments, Jc is identified as the current density J at 
which the electric field E exceeds a threshold value. Usually pinning prop- 
erties of the vortex lattice are uniquely characterized by the Jc value and the 
corresponding E — J relation. However, in the peak region of 2H — NbSe2, 
Jc is found to increase significantly with the rate at which J is increased!^. 
This amounts to the Jc vs T curve being much sharper when a large value 
of dJ/dt is used|^]. Moreover, E — J curves themselves exhibit significant 
hysteresis in upward and downward current cycles Q, |5|, y]. 

The history dependence of Jc under different quasi-static field excur- 



sions was recently explained by a phenomenological model 1 12]. This model 
supposes that there exists a set of metastable vortex configurations, each 
corresponding to a different value of Jc- At a given field and temperature, 
different configurations differ only in the extent of lattice order. Jc can thus 
be viewed as a macroscopic representation of a particular metastable vor- 
tex configuration. The most important assumption of this model is that 
there exists a stable vortex configuration with Jc = Jc^ which is unique for 
a given field and temperature. Depending on the particular magnetic his- 
tory, a metastable configuration with Jc either smaller or larger than J^* 
is produced. The extent of history dependence is assumed to be governed 
by the parameter Br, which is a macroscopic measure of the metastability. 
Evolution of Jc when the field changes from B to B + 6B is postulated to 
be governed by 



UB + 5B) = UB) + {\5B\/Br)[r,\B) - J,(S)] (1) 

Eq. 1 provides that a metastable configuration with Jc 7^ J^* would evolve 
into the stable state when the field is repeatedly cycled by a small am- 
plitude. Recently, the existence of such a stable state has been confirmed 



experiment ally 1 13]. Physically, we can imagine that the process of repeated 
field cycling pumps the vortices in and out of the superconductor. In the 
absence of adequate thermal energy field cycling allows the vortex lattice 
to explore the configuration space and reorganize or anneal into the most 
stable vortex configuration. Thus the process of annealing is aided by the 
motion of vortices caused only by an external driving force such as a field 
change or a transport current larger than Jc^ 14]. 



We believe that the dynamic phenomenon mentioned above and the his- 
tory dependence in critical currents are two different manifestations of the 
same phenomenon of metastability in the vortex state and therefore should 
be described by a common phenomenology. In this paper, we present a gen- 
eralization of the phenomenological model described above, to account for (i) 
the history dependent dynamic response of the vortex lattice in the peak re- 
gion (c.f Fig. 4 of Ref. J^); (ii) the hysteresis vnE—J curves ]|^, ^, P; (iii) the 
dependence of Jc on the current sweep rate used m E — J measurements P] 



and (iv) the so called negative dynamic creep phenomenon] 10 1. 

We postulate that, the time evolution of a metastable vortex configura- 
tion in terms of Jc is described by the equation 

dJc/dt = (Jf - Jc)/t (2) 

where t is the time variable and r is the time constant for the annealing 
process. We argue that the time constant r is inversely related to the velocity 



V of vortex motion, i.e., r ~ \v\~^ ~ ]-^]~^114, [l^. In other words, a larger 
velocity of vortices facilitates faster rearrangement of vortices into the stable 
configuration (annealing). The absolute value of E (or v) implies that the 
process of rearrangement is independent of the direction of vortex motion. 
Let us first examine the vortex rearrangement driven by a quasi-static 
field sweep, in the absence of transport current. If the field sweep rate is 
dB/dt, E ~ dB/dt and the time constant r (~ ]^I^^) can be expressed as 
TQR/\dB/dt\ where Ris a characteristic value of the sweep rate below which 
the field change can be considered quasi-static. The time constant tq{B,T) 
can be a function of the field B and temperature T. Substituting the above 
expression for r in Eq. 2, we get 



dJJdB = ±(jf - Jc)/Br (3) 

where B^. = tqR with upper (lower) sign for the field increasing (decreas- 
ing) case. Eq. 3 is the differential form of Eq. 1 describing the magnetic 



history dependent Jc [ 12 1. Therefore, we conclude that the parameter B^ 
which is a measure of the history dependence in Jc, is directly related to the 
time constant tq which, as would be shown below, governs the dynamics of 
metastable vortex states subjected to an external driving force. 

We now consider the dependence of Jc on magnetic field sweep rate for 
a fixed applied field. When the external field B(t) increases (or decreases) 
linearly in time, we can write t = ±[B{t) — Bi]/B. The upper (lower) sign 
is again applicable for increasing (decreasing) field. Bi is some initial field 
where the field sweep started and B (positive) is the magnitude of the sweep 
rate dB/dt. Substituting for t in Eq. 2, we get 

dJc/dB = T[{B{t) - Bi)/B\j^' - Jc)/t (4). 

In the increasing field cycle (upper sign) {B{t) — Bi) is positive and Jc < Jc^ 
in the peak region[^, 13|. Therefore dJc/dB < (negative dynamic creep 



as observed in experiments |1C]. On the other hand, in the decreasing field 
cycle (lower sign) dJc/dB > because Jc > J** (see Ref. Q) and {B{t)—Bi) 



is negative. In other words, the negative dynamic creep observed in the 



increasing field cycle [ 10 1 is a manifestation of the finite time delay needed 
for the vortex lattice to become more disordered. 

We now focus on the transport current driven annealing process governed 
by Eq. 2 where the time constant r ~ \E{J)\^^. For E{J), we choose the 
simple form E{J) ~ {J/Jc — 1) for J > Jc and thus 

r(J) = To/(J/Jc - 1) (5). 

For J < Jc, E{J) = and therefore r — > oo. The specific form chosen for 
E{J) is however not critical for the results we are going to present. Any 
other form which makes r —> oo as J ^ Jc would equally serve to illustrate 
the results in this paper. 

Eq. 2 and Eq. 5 suggest that irrespective of whether (i) Jc < Jc^ or 
(ii) Jc > Jc^, a transport current J > Jc causes Jc to increase or decrease 
respectively towards the stable value. In other words, while J > Jc, Jc and 
thereby r change with time, resulting in the time dependence of the voltage 
response through the relation E{J) ~ {J / Jc{t) — 1), even when the transport 
current J is kept fixed. 



Combining Eq. 2 and Eq. 5 and substituting the dimensionless quan- 
tities jc{t) = Jc{t)/J^^, j = J/Jc^ and t — > t/TQ, a formal solution for jc{t) 
can be written as 

Ut) = 1 + [je(to) - 1] exp[- / dt{j/Ut) - 1)] (6) 

where jc(io) is the normalized critical current of the vortex configuration in 
which the system is initially prepared at time to- Here we have two cases of 
interest viz., (i) jdto) < 1 and (ii) jc{to) > 1 which are henceforth referred 
to as superheated and supercooled vortex configurations respectively! 13 1. The 



superheated configuration can be achieved by entering the peak region either 
by increasing the field at constant temperature 1 12, 13] or by increasing tem- 



perature at a constant field. On the other hand, the supercooled state is 
achieved by either decreasing the field from above the normal state bound- 



ary while keeping the temperature constant [^, |l^, 13 1 or by cooling the 
superconductor in a field 0, ^. 

We solve Eq. 6 iteratively to obtain the time evolution of jc(t) in response 
to a step increase in the transport current density j at t = tg- In Fig- 
1, we present the time dependence of jc{t) and E{t) ~ {j/jc{t) — 1) for 
the superheated vortex configuration, i.e, Jdto) < 1- The voltage response 
decreases with time as observed experimentally (c.f. Fig. 4(b) of Ref. |^] 
and Fig. 2 of Ref. ||9|). In the inset of Fig. 1(b), we show E{t) vs ln{t) at 
different driving currents, which is in excellent qualitative agreement with 
the experimental results (c.f. Ref. Q). In Fig. 2, we present jc{t) and E{t) 
obtained by starting from a supercooled vortex configuration (i.e., Jdto) > 1); 
which is again in good qualitative agreement with experiments(c. f. Fig. 
4(b) of Ref. |]). 

Let us now consider the situation in Fig. 1, but with the driving current 
j switched off to zero during the interval ti < t < t2- At time ti, the vortex 
lattice is frozen into a metastable configuration corresponding to jc = jci^i) 
as the time constant r ^ oo for j < jc- When the current is switched on 
again at t = ^2, jc begins to evolve from the value jc{ti) according to Eq. 
6 as shown in Fig. 3(a). Similarly, one can consider the case of reducing 
the current j by a small amount over a time interval as shown in Fig. 3(b). 
This merely slows down the process of evolution of Jc towards a stable value 
J**. The calculated voltage responses shown in Figs. 3(a) and 3(b) compare 
very well with the results of Ref. |^ . 

The above discussion provides a framework for the annealing process 



mediated by a field change [12, T^. When the critical state is set up, shielding 



currents (just like eddy currents) are set up at a high current density (>> 
Jc) ||l6[| which decay to the value Jc with a characteristic time constant T ~ 
Hoo^/Apff governed by the flux flow resistance pjf = pnH/Hc2, where p„ 
is the normal state resistivity, H and Hc2 are the external field and the 
upper critical field respectively, ^o = ^vr x 10~^ W/A.m and a is the sample 
dimension. This is equivalent to a pulse of transport current Jmean > Jc of 
duration T which causes a partial annealing of the metastable vortex state. 

This model can also be used for calculating the E — j relations. We 
consider that the transport current density j (> jc) is increased in small 
steps of Aj each of duration At. We calculate the jc and the corresponding 
E{j) values at the end of each current step. In Fig. 4 we plot the E — j 
curves obtained with two different values of the current sweep rate dj /dt 
(~ Aj'/At). Clearly, the critical current density inferred from these curves 
would be smaller when dj /dt is larger, just as observed in experimental 
studies(c.f. Fig. 1 of Ref. Q). Further, the shape of E — j curves differ 
dramatically with the current sweep rate. 

We also studied the E — j curves in both upward and downward current 
cycles with a fixed dj /dt. In Fig. 5(a) we consider the case where the vortex 
lattice is iniatially prepared in a metastable configuration with jdto) < 1 
{superheated state) and in Fig. 4(b), for the case jc(to) > 1 {supercooled 
state). The calculation clearly brings out the hysteretic E — J behavior 
usually seen in the peak region and the results compare remarkably well 
with the experiments IQ, ^, |6|. We further note that, after current cycling, 
jc (Jc) always settles to the value 1 {Jc = Jf) as shown in Ref. [§. 

The results presented here are contrary to the notion that the E — J 
relation and the Jc value inferred from it, uniquely characterize the pinning 
properties of the underlying vortex structure in all circumstances. We ar- 
gued instead that while the transport current J is higher than the threshold 
value Jc, the Jc value itself changes irreversibly, leading to the hysteresis in 
the E — J relations. 

We caution here that the simple model presented here explains a variety 
of dynamic response of the vortex lattice to a transport current above the 
macroscopic threshold value Jc- However, it can not account for the ac 
dynamic response for current values below Jcfl^l' ^ qualitative mechanism 
for the ac dynamic response was recently presented by Paltiel et al[|l4|] by 
considering the disordered vortex phase, injected into the superconducting 
sample across the surface imperfections, coexisting with the ordered phase 
in the interior of the sample. However, an explicit connection between this 
mechanism and the model developed here is not yet available. 



In conclusion, we have presented a macroscopic model which accounts 
for the dynamics of the vortex lattice in the peak effect region. We ex- 
plained the history dependent dynamic response of the metastable vortex 
states subjected to a transport current and the hysteretic voltage-current 
characteristics observed in the peak region of weakly pinned superconduc- 
tors such as 2H — NbSe2- We are also able to explain the strong variation 
in voltage current characteristics with the current sweep rate. Further, the 
so called phenomenon of "negative dynamic creep", viz., the anomalous de- 
pendence of Jc on the magnetic field sweep rate in the peak region is also 
explained. 

The author thanks Prof. S. Bhattacharya, Dr. K. V. Bhagwat, Prof. A. 
K. Grover, Prof. S. Ramakrishnan and Dr. V. C. Sahni for discussions and 
critical reading of the manuscript. 
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FIGURE CAPTIONS 

Fig. 1: (a) Time evolution of jdt) in response to a step change in 
the transport current j (dotted line). Initial vortex state corresponds 
to a superheated vortex configuration with jdto = 0) = 0.5 (< 1). 
The corresponding voltage response is shown in (b). In the inset of 
(b), this response is plotted as E vs ln{t) for different drive currents 
(compare with Fig. 2 of Ref . M ) . 

Fig. 2: Time dependence of jc{t) and E{t) in response to a step 
current pulse (dotted line). The initial vortex state corresponds to a 
supercooled vortex configuration with jc(io = 0) = 1.5 (> 1. 

Fig. 3: (a) Time dependent response of the superheated state consid- 
ered in Fig. 1, but with the drive current j switched off during the 
interval ti < t < ^2- Note that E{t = ^2) = E{t = ti) which is referred 
to as the "long term memory" effect (see text), (b) same as (a) but the 
current j (dotted line) is only lowered by a small amount (j > Jdt))- 
Note that the time dependence of E{t) is significantly slowed down. 

Fig. 4: Calculated E — j relations for two different values of current 
sweep rate dj/dt as indicated, j* = 0.5 is the critical current of the 
vortex state in which the superconductor is initially prepared. 

Fig. 5: Calculated E — j relation for upward and downward current 
cycles as indicated by the arrows. In (a) vortex state is initially pre- 
pared with a critical current density j* < 1 and in (b) initial state 



corresponds to j* > 1. The dj /dt is kept fixed at 0.1. In both cases 
-E(j) becomes zero at j = 1 {Jc = J^*) i'^ the downward current cycle. 
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